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I. INTRODUCTION 



In the limit of infinite Coulomb repulsion (U = oo ) 
one can transform the Periodic Anderson Model (PAM) 
with two channel (SU(2)) r hy employing the Hubbard 
Xj ta b=\j,a) (j,b\ operators,!!! where the set {\j, a}} is an 
orthonormal basis in the space of interest. Projecting 
out the components with more than one electron from 
any local state at site j one obtains 



H = J2 E K J^c Ka + E fd* X . 
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The first term is the Hamiltonian of the conduction 
electrons (c-electrons) . The second term describes in- 
dependent localized electrons (/-electrons), and a sim- 
ple index j is used to indicate the sites. The last 
term is the hybridization Hamiltonian giving the inter- 
action between the c-electrons and the /-electrons with 



= (l/\/A7)V tT (k) exp (ik.Rj), where Rj is the po- 
sition of site j and N s the number of sites. Notice that 
this interaction conserves the spin component a. 

The X operators do not satisfy the usual commutation 
relations. One has then to use the product rules: 



been found that it is not usually satisfied when the n a 
are calculated with approximate cumulant Green's func- 
tions (GF)b An approximation displaying this behavior 
is the "Chain Approxjiaation'' (CHA), which was first 
employed by HewsonJj'El and is the most general cumu- 
lant expansion with onkj-jsecond order cumulants, as well 
as being $— derivable .Qu In a previous worktl we have 
stated a conjecture that gives a systematic way of achiev- 
ing completeness by adding a set of diagrams to an ar- 
bitrary family. This result was verified for the CHA, 
and the resulting set of diagrams is not derivable any 
more. 

In the present work we recover completeness by a pro- 
cedure similar to that employed by ColemanB in the slave 
boson approach. This last method predicts unphysical 
second order phase transitions both for all parameters 
at intermediate temperatures T and in theiKondo region 
(fx >> Ej\j a = Ef) at low T. ColemanlI2l has observed 
that these effects are artifacts of the theory, and the ad- 
vantage of the present treatment is that those spurious 
phase transitions do not occur. Our method gives results 
very close to those obtained by the slave boson in the 
Kondo limit at low temperatures, while it recovers those 
of the CHA at high T for all parameters. 



II. X-BOSON CUMULANT METHOD 



j,adt 



(2) 



and the diagrammatic methods based on Wick's theo- 
rem are therefore not applicable. We shall use instead 
a cumulant expansion that was originally employed by 
HubbardB to study phis model, and that was latter ex- 
tended to the PAMBH 

There are four local states at each site j of the lat- 
tice: the vacuum state | j, 0), the two states | j,a) of 
one electron with spin component a, and the state | j, 2) 
with two local electrons. When U — > oo the state | j, 2) 
is empty, and we have used the Hubbard operators to 
project it out from the space of local states at site j. In 



this space, the identity 
the relation: 



Ij at site j should then satisfy 



X 
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where a is the spin component opposite to a, and the 
three X 3 ^ aa are the projectors into | j,a). The occu- 
pation numbers n^ a =< Xj >aa > can be calculated from 
appropriate Green's functions (GF), and assuming trans- 
lational invariance we can write rij^ a = n a (independent 
of j), so that from Eq. ([|) it follows that 

n + n a + n w = 1. (4) 

We shall call this relation "completeness", and it has 



The—present work is a modification of a preliminary 
versiorcll that was inspired in the mean-field-approxima- 
tion of Coleman's "slave boson" method.00 He writes 
the Hubbard X operators as a product of ordinary bosons 
and fermions: Xj >o0 — * b+bj, X jt0a — ► b+f JjCr , Xj,ao — ► 
fj^ a bj, and uses the equivalent of our Eq. (Q) to avoid 
states with more than one electron at each site j. In the 
spirit of the mean field approximation bf b^ >= r, 
and the method of Lagrangian multipliers is then used to 
find the "best" Hamiltonian that satisfies Eq. (||). The 
problem is then reduced to an uncorrelated Anderson lat- 
tice with renormalized hybridization V — > rV and / level 

e / e f + x - 

Our method consists in adding the product of each Eq. 
(||) times a Lagrange multiplier Aj to Eq. (]l|), and the 
new Hamiltonian generates the functional that we shall 
minimize employing Lagrange's method. Instead of the 
parameter r we introduce 



R 



E < Xj. 00 >- 



-< bp 3 >, 



(5) 



and we call the method "X-boson" because the Hub- 
bard operator Xj. oa has a "Bose-like" character.!!! Con- 
sidering Eqs. (^,^), and employing a constant hybridiza- 
tion V, as well as site independent local energies Et : j t(7 = 
Ef t(T and Lagrange parameters Aj = A, we obtain a new 
Hamiltonian with the same form of Eq. (m 
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h = J2 E ^ 4, 

k,<r 
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but with a renormalized localized energy 



A. 



(6) 



(7) 



This procedure allows for an independent variation of 
R when the Free Energy is minimized, even though the 
completeness relation 



R = 1 - < X ° 



> 



(8) 



must be simultaneously satisfied. 

In the slave-boson approach a one-body Hamiltonian 
is obtained at this stage, and Eq. (Q) is then necessary to 
avoid the occupation of more than one localized electron 
per site, while the conservation of probability in the space 
of the local states is automatically satisfied, because nor- 
mal fermion operators are employed in the transformed 
Hamiltonian. The Eq. (Q) on the other hand employs 
X-operators, that force the local states to be singly oc- 
cupied at most, but Eq. (Q) (completeness) must be im- 
posed here because it is not automatically satisfied when 
approximate GF are employed to calculate the n a . In the 
present work we shalLemploy the GF of the Chain ap- 
proximation (CHA), Era because they give a fair descrip- 
tion of the system in spite of their simplicity. 

The present treatment employs the Grand Canonical 
Ensemble of electrons, and instead of Eq. (fy) we shall use 



Ti = H - [1 ^ ^ C^Ck.a + ^2 V aX haa 
k,a" ja 



(9) 



where v a = 0, 1 is the number of electrons in state 
a >. It is then convenient to define 



and 



£j,a — Efj^a — [W a 



£kcr — Ek a — A* 



(10) 



(11) 



because Ef t j ia and appear only in that form in all 
the calculations. The exact and unperturbed averages of 
the operator A are denoted in what follows by < A >n 
and < A > respectively. 



III. THE CHAIN APPROXIMATION GREEN'S 
FUNCTIONS 

The CHA gives simple but useful approximate propa- 
gators, obtained in the cumulant expansion by taking the 



infinite sum of all the diagrams that contain ionic vertices 
with only two lines. The laborious calculation of the gen- 
eral treatment!!! is rather simplified in this case, and we 
shall give a brief description of the technique, particularly 
when only the imaginary time is Fourier transformed, 
because this intermediate situation is not discussed in 
reference |l| and it is necessary to calculate the impurity 
problem. The only X a and X\ operators of the Fermi 
type that appear in the calculation have a = (0, a), and 
the f-electron GF in real space and imaginary time are 



G ff (j,a,T;jr,af,T/) 



n(i2) 



where Xj t0 ,(r) = exp (rTi) Xj iOC exp (— rTi) corresponds 
to the Heisenberg representation and the subindex + 
indicates that the operators inside the parenthesis are 
taken in the order of increasing r to the left, with a 
change of sign when the two Fermi-type operators have 
to be exchanged to obtain this ordering. El In a similar 
way one defines the c-electron GF G cc (k, a, r; k/, a/, r/) 
as well as the mixed GF G^ c (j, a, r; k/, af, ti) and 
G cf (k, a, r;jl, al, Tl). 

The boundary conditions of these GF with respact to 
r makes possible to expand them in Fourier seriesEI em- 
ploying the Matsubara frequencies u) n — (2n + l)i7r//3, 
where n is any integer. Because of the invariancc of Ti. 
against r translations we have frequency conservation,!!! 
so that Gff(J, a, iu) n ;jl, a/, iu) n l) = unless u) n +ui n l = 0, 
and one can then write 

G ff {j, a, iun'Jf, at, iw n f ) = G jl-j/ a >( z n) A(w„ + w n t) 

(13) 



and similar relations for the G CC ,G-> C , and G c * . Here we 
have used z n = iui n , and we shall keep this notation in 



what follows, as well as employing A (7) = 5 7i o (& modi- 
fied notation for Kronecker's delta). The "bare GF" (i.e. 
with all V,,k,cr = 0) take a fairly simple form because 
many of the relevant operators become statistically inde- 
pendent, and we then have 

G ff (j, a,z n ; j/, a/, z n i) =*> 

G °f,a( Z n) A(w„ +UJn') Sjj, 6 

a, aft 

(14) 

G cc (k, a, z n ;k/, a/, z n i) 

Gc.koX 2 ") A(w„ +uj n i) <!>k,k/ S a ^ h (15) 
G fc =*> G cf =>■ 0. (16) 

The <5 ajQ / in the bare GF follows from the commutation 
of Ti with the z component of the spin. 

The diagrams in real space that contribute to the CHA 
are schematized in Fig (Q). The meaning of the symbols 
in the cumulant diagrams is the following: 

a) the "vertex" % = G° f 0a (z n ) = -D 0<7 /(z n - £/ jCT ) 
is the / bare cumulant GF , where 

Doa —< X oa > + < X aa > 

b) the "vertex" Q = G° c ^{z n ) = -l/(z n - e k(T ) 
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is the c bare cumulant GF 

c) the lines (edges) determine an open loop with a def- 
inite direction. When the line points to the f-vertex is 



< — =Vj,k,<r, while < — = ^,*k a when it points to the con- 
duction vertex __ 

d) The cumbersome sign and symmetry factorsEIn are 
rather simple in the CHA: it is only necessary to mul- 
tiply the G^ c and the G c f into a minus sign. 

e) As both the uj and the a are conserved along the open 
loop, it only remains to sum over all the internal j and 
k. 




ki 



(18) 



which is the local GF at site j times \V&\ when the hy- 
bridization is purely local, i.e. when V a (k) — V a . The 
Eq. ([j"?]) is then a geometric series that is easily summed: 



GpavM = S °>,°hj 1 _ G oj Zn)M<r{Zn) (19) 

_A)CT (20) 



£/> + D 0a M a {z n ) ' 



Employing the same technique, we obtain for the CHA 
the G cc , that gives the scattering by the local impurity 
of a conduction electron k,cr into k/,e>7: 



ja ka 

FIG. 1. The Green's function diagrams in the CHA. The 
filled circles (vertices) corresponds to the /-electron cumu- 
lants and the empty ones to those of the c-electrons. The 
lines (edges) joining two vertices represent the perturba- 
tion (hybridization) (cf. Section |tll|) a) Diagrams for the 
/-electron GF in the CHA. The CHA diagram for the 
/-electron Gjf a .j faf (iLu n ) is represented by the filled square 
to the left, b) Diagrams for the c electron GF in the CHA. 
The G£5 - /; k, CT (- z ii) i s represented by the square symbol to the 
left, c) Diagrams for the f-c electron GF in the CHA. The 
Gj.£ CT (zn) is represented by the square symbol to the left. 

We then obtain (with a — (Oct)): 

Gji a r,ja( Z n) = &al ,aG°f a {z n ) 

+G° fa (z n ) 2J Vj'McrG° Mia (z n )V* kicr G f a (z n ) 

ki 

+G°f a (z n ) Vj/facrG" k 2<J (z n ) 22 Vj 1 te<jG°f^ a (z n ) 

k 2 ji 

x V JlM aG° cMa (z n )v; Ma Gl a (z n ) + (17) 



the sums over the internal j r reduces to a single term 
with jl = j, and we introduce 



Jk,k/ 



(21) 



Gkral:kcr{ z n) — &al,a { G° k(T (z n ) (5k 

+G° k , ff (z„)V-;(k/)Gj4 Q (z„)xy ff (k)G° kCT (z„)}, (2 . 
as well as the mixed GF (with a' = (Oct')) 

where we have already included the minus sign (discussed 
in rule d) above) into 

= -G f / a . ja {z n )V jM G° cM {z n ) (23) 



(24) 



2. The lattice case 

The case of the GF in reciprocal space and imaginar 
frequencies has been discussed in detail in reference,! 
and in the CHA one follows the same prescriptions given 
above, but replacing the sum over internal j and k by a 
conservation of k along the whole chain, so that we have 



1. The impurity case 
When there is a single impurity at a given j, each of where 



G ff (k/,{0af),z n f; k, (0a),z n ) 

= SuMl $a,al A(W„ + W„/)G^(z„), (25) 

G cc (k/, at, z n l; k, cr, z n ) 

= <J k ,k/ 6„ ttrt A(w„ + uj n f)Gi c a (z n ), (26) 
G /c (k/, (Oct/), z n l; k, cr, z n ) 

= (5k.k/ <J CT)0 -/ A(w„ + w„/)G^(z„), (27) 



GH(Zn) = 



(z n - £f,a) {Zn - £k CT ) - |K(k)| 2 Ar' 



(28) 
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G^ c a {z n ) 



~ On - £/,« 



and 



{Zn ~ £f,a) (Zn ~ Ska) ~ \V a {k)\ 2 D a ' ^gs 

- D a V a (k) 



(z n - e />CT ) (z„ - e k(T ) - |Vv(k)| 2 D CT 



(30) 



5. The Slave Boson GF vs the GF of the CHA 

The Slave Boson GF in the mean field approximation 
(cf . references |l2"|-|l4]) are just the GF of the uncorrelated 
problem (U — 0) but with a renormalized hybridiza- 
tion V — > V = rV plus a condition that forces the lo- 
cal electron to an occupation less or equal to one. The 
GF of the CHA given above are formally very close to 
the uncorrelated ones, but they can not be reduced to 
them by any change of scale (except for D a = 1, when 
we recover the slave-boson Green's functions if we put 

V -► rV = V in Eqs. dHJ^) and Eqs. (|ffcj)). The 
obvious change ^/Dq^V — > V, leaves an extra factor D$ a 
in the G f J at . ja (z n ) and G(£(z n ), as well as a VAkr in 

both G^ a (z n ) and G^ r (z n ), and these factors are re- 
sponsible for the correlation in this approximation, aad 
lead to essential differences with the uncorrelated case.tS 
In particular, they force the total occupation rif of the 
f electron to n/ < 1, while in the uncorrelated case the 
relation rif < 2 is satisfied. In the slave boson method 
the imposed condition rif < 1 is fulfilled by a shift in 
the local energy e^ a — ► £y )Cr = e/ -cr + A and a reduction 
of the hybridization to V —>■ V = rV . From an opera- 
tional point of view, a shift in s / j(T might not be sufficient 
to force rif < 1 because the hybridization extends the f 
spectral density to the whole conduction band, and re- 
ducing V helps to satisfy this condition. By increasing 
the temperature T or the chemical potential \x the value 

V = is presently reached, leading to a decoupling of 
the two type of electrons that can be interpreted, from a 
more general point of view, as a change of phase related 
to a symmetry breaking of the mean field Hamiltonian. 
Although the condition that forces completeness in the 
CHA is identical to that employed in the slave boson 
method to force rif < 1, it has a rather different origin, 
being only a consequence of using a-.reduced set of dia- 
grams in the perturbative expansion,!!! and the departures 
from completeness are usually very moderate. In the for- 
malism described in the present work, it is this essential 
difference between the two methods that eliminates the 
spurious phase transitions appearing in the slave boson 
method. 



IV. THE SINGLE IMPURITY PROBLEM 



In the X-boson approach D a = R + nt a must be cal- 
culated self-consistently through minimization with re- 
spect to the parameter R of an adequate thermodynamic 
potential. When the total number of electrons N t , the 
temperature T and the volume V s are kept constant one 
should minimize the Helmholtz free energy, but the same 
minimum is obtained by employing the thermodynamic 
potential il = — fcsTln(Q), (where Q is the Grand Par- 
tition Function) and keeping T, V s , and the chemical 
potential [i constant (this result is easily obtained by em- 
ploying standard thermodynamic techniques). 

A convenient way of calculatiaa-XLis to employ the £ 
parameter integration mcthodJi3'E3'EZI This method in- 
troduces a £ dependent Hamiltonian if (£) = ff + £ifi 
through a coupling constant £ (with < £ < 1), where 
ifi is the hybridization in our case. For each £ there is 
an associated thermodynamic potential which sat- 
isfies: 



dfl 



--<Hi(0 >C, 



(31) 



where < A >^ is the ensemble average of an operator A 
for a system with Hamiltonian if (£) and the given values 
of /i, T, and V s . Integrating this equation gives 



= n o + f de<ffi(0 ><;, 

Jo 



(32) 



where Q a is the thermodynamic potential of the system 
with £ = 0. This value of £ corresponds to a system 
without hybridization, and one obtains (in the absence 
of magnetic field = £k an d £/er = £/) 



n o = -~^ln [1 + exp(-^ k )] 



In [1 + 2 exp(-/3g»] + A(R - 1), 



(33) 



and to calculate f2 in Eq.(p2j) we use 



< H x >£= 2Re 



E 

kcr 



(34) 



The average yc ko .Ao CT y is obtained from the ana- 
lytical continuation of the Matsubara C?/ k a { z ni £) —> 
f c 

Gj. ko .(z,£) into the complex upper and lower semi- 
planes, where G^ kcr (z„,£) is the GF in Eq. (23) but for 
Vj,k,a— ^V^k.o-- One then finds 



oo 

— / duj npluj) 



(35) 
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where n F (x) = 1/ [1 + exp((3x)] is the Fermi-Dirac dis- 
tribution. From Eqs. (|l8| , p4| ,35) we then obtain 



(Hi)t = - ^ J du n F (uj) 



ZD 0a M a (u+) 



u+-e f +D 0a ?M <T (w+) 



(36) 



where lu + = lu + iO. This equation has an interesting 
scaling property: it is equal to the corresponding expres- 
sion of the uncorrelated system for the scaled parame- 
ter V^kjcr = y/Doo-Vjfcv (it is enough to remember that 
by replacing D$ a = 1 in the GF of the CHA one ob- 
tains the corresponding GF of the uncorrelated system). 
Rather than performing the £ and ui integrations, we shall 
use the value of the Vl u for the uncorrelated system with 
Vj,k,a = yfDoaVjfco an d employ Eq. ( |32| ) to calculate 
Jo d£, < £TJ»(f) >£= n u - fi« where 



fi" — — 

+A(R-r 



In [1 + exp(-^e k )] + In [1 + exp(-/3e»] 

(37) 



L k 



where 



fin 



In 



l + 2exp(-/3g» 
(1 + exp(-/?£ / )) 2 



(42) 



and the eigenvalues iai^ a of the TL U are just given by the 
poles of the GF in the' CHA (Eq. (ph). 



It is interesting that all the correlation effects on the 
thermodynamic potential appear in the f2 , that corre- 
sponds to the unperturbed system. The effect of the 
perturbation in the present approximation is to redis- 
tribute the energy of the quasi-particles in the same way 
as in an uncorrelated system with hybridization constant 
Vj,k,<? = y/DotrVjfce, and one could then expect Fermi 
liquid behavior in the CHA (see a formal discussion in 
the Appendix). 

The parameter A is obtained minimizing Q with re- 
spect to Rt3 (at constant fi as discussed at the beginning 
of this Section). To simplify the calculations we shall 
consider a conduction band with a constant density of 
states and width W = 2D, and an hybridization con- 
stant Vo-(k) = V. We then obtain 



dR 



~dR 



+ A = 0. 



(43) 



is the Vt u for Vj t k 1<7 — 0. In our case the unperturbed 
Hamiltonian for the lattice problem is 



k.tr 



f}. Jj., + N s K{R-l) 



3,v 



(38) 



and in the impurity case the sum over sites reduces to the 
impurity site and N S A(R - 1) -> A(R - 1). This Hamil- 
tonian can be easily diagonalized, and the corresponding 
H u can be written 

n u = J2 a L a ^ + A (R - i), (39) 



where a\ a (cti^) are the creation (destruction) operators 
of the composite particles of energies cjj i<7 (there are N s + 
1 states for each spin a). The calculation of 



n u = -j ln [1 + ex P (-/3 Ui,a)] + A(i? - 1), 



(40) 



is straightforward, and employing Eqs. ( |32| , |3"3| , [37 40) we 
find 



fi = fi + — Y ln [1 + exp(-/3 uj ha )} + A(R - 1) 



(41) 



The poles of the impurity GF satisfy: 



£ f 



V 2 Dg 

2D 



In 



D 



D 



(44) 



and calculating the R derivative of u>i(k, a) we obtain 



A 



-V 



In 



a + D\ D 2 )n F (u^) 



2D ^ \u) i>a -Dj {u? a -D 2 + V 2 D, 



j (45) 



In the absence of magnetic field, i.e. when the Wj )0 - = 
is independent of a, we can write this equation as: 



A ~ V r A I M f uj + D 

A = / dio Pfiw) in — 

D /_„ m ' \w-D 



' — OO 

{lu 2 -D 2 ) n F (m) 
( w 2 _ D 2 + V 2 D a ) ' 



(46) 



where Pf(u>) =Y2i^( LJ ~ ^i) ^ s the density of / states, 
and the chemical potential /j, of the electrons was included 
into the single particle energies (cf. Eqs. ( ]To| , pl] )). The 
Kondo temperature 



-K 



D exp I -^4 



(47) 



is defined as the T that makes the slave boson parameter 
r to vanish, but another useful definition was given by 
Bcrnhard and LacroixtLa by taking Tk as the crossover 
temperature, determined by the maximum of the deriva- 
tive of (el fia) with respect to T. 
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V. IMPURITY RESULTS 



The results presented in this section correspond to 
a Kondo regime with the following parameters: Ef = 
-0.15; W = 2.0, V = 0.3. 

Figure |^ shows the evolution with T of the parameter 
R, that measures the hole occupation, together with the 
corresponding r 2 of the usual mean field slave boson, as 
well as the two occupation number n j . The figure shows 
that r 2 — ► at a finite T (Kondo temperature Tk) 
while R remains positive, avoiding the spurious phase 
transition. The n/ — > 1 at Tk for the slave boson, while 
it remains lower than 1 for the X-boson at all T. 

At high temperatures the system is well described by 
localized moments coupled to the conduction electrons 
through a Coqblin-Schriefer type Hamiltonian, and even 
at high temperatures there should not be a complete 
decoupling nor a second order phase transition, as sug- 
gested by the mean field slave-boson theory.Ea As the 
GF of the CHA retain the correlation properties of the 
system, the X-boson approach shows no decouplin g nor 
spurious phase transitions, as discussed in Section III S . 
One can see in Figure || that the two approaches give 
similar results at low temperatures. 



1.0 



Slave boson / 
X-boson / n ' 




0.0 



10' 



10 



10" 



10 



FIG. 2. Occupation numbers n f , and parameters r and R 
as a function of T for both the slave boson and the X-boson 
methods. The figures correspond to the following parameters: 
E f = -0.15; W = 2.0; V = 0.3; fx = 0.0. 

In Figure [| we represent the parameters A and A as a 
function of temperature. We observe that the slave boson 
A breaks down at the Kondo temperature Tk whereas the 
X-boson A goes continuously to zero. In the high tem- 
perature limit the results obtained with the usual CHA 
are recovered by the X-boson. In Figure || we show the 
evolution of A, r 2 , and n/, as a function of the chemical 
potential fi in the usual slave boson approach: the for- 
malism breaks down at a value /in, where n/ — > 1. In 
Figure [| we show the evolution of the parameters A, R, 
and n/ as a function of the chemical potential [i for the 
X-boson approach. The results of the two approaches are 



similar for /i < /in., i.e. before the breakdown of the slave 
boson method (compare Fig. | and Fig. |). The param- 
eter A reaches a maximum in the Kondo region and goes 
to zero when fj, » Ef, where nj -> 1, and in this limit 
we recover the usual CHA. 




FIG. 3. Evolution with temperature T of the renormalizing 
A and A, for the same model parameters used in Figure 0. 




FIG. 4. Slave boson: A, r 2 , and n/ vs. /i for the same 
model parameters used in Figure ^| and T = 0.001. 




FIG. 5. X-boson: A, R, and rif vs. fi for the same model 
parameters used in Figure H. 
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The values of the f and c electron density of states 
Pf(p), and pcO-O at u> — p are shown as a function of p 
for the slave boson method in Figure || and for the X- 
boson approach in Figure [?]. The slave boson plot breaks 
down at p — po ( where rif — > 1) while for the X-boson 
the density of states Pf{p) is maximum in the Kondo 
region and goes to zero when p >> Ef. In both cases 
one observes the transfer of conduction electrons to the 
/ band giving rise to the Kondo resonance. 
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FIG. 6. Slave boson: Pf(u), Pc{p) vs. for the same model 
parameters used in Figure 




FIG. 8. Density of states of Pf(oj), p c (w) vs. ui in a typical 
Kondo situation, in the two approaches, for the same model 
parameters used in Figure H and fi = 0.0. 



VI. THE LATTICE PROBLEM 



For the lattice we follow the same technique employed 
with the impurity: the parameter A is obtained by min- 
imizing the thermodynamic potential SI = — fcsTln(Q) 
with respect to R and D a = R + nj a is calculated self- 
consistently. The f2 is again obtained from Eq. (^) with 
adequate values of fin and < 
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FIG. 7. X-boson: Pf{p), Pc(p) vs. p for the same model 
parameters used in Figure 0. 



In Figure || we present the density of states Pf{uj) and 
p c (w) vs. energy u> in a typical Kondo situation, both 
for the slave boson and X-boson treatments. The Kondo 
resonance appears in the two cases with a similar shape, 
but it is less pronounced in the X-boson treatment. 



f> = --]Thi[l + exp(-/?£ k )] 



In [1 + 2exp(-/3S»] + N S A(R - 1), (48) 



(Hi) 



k.o- k 



£ \v4V\ 2 D a 



(oj+ - e ta ) (oj+ - e kCT ) - e |K(k)| 2 A, U 



49) 



where c<j + = uj + iO. As in the impurity case, this ex- 
pression coincides with that obtained for an uncorrelated 
Hamiltonian, which in this case is Eq. ( ^8|) with a scaled 
hybridization V a (k) = v / A^V £T (k). The corresponding- 
eigenvalues can be calculated analytically in this case, 
because for each spin component a the Hamiltonian is 
reduced into N s matrices 2x2. They are the poles of the 
lattice GFs, replacing the N s + 1 eigenvalues of the 
impurity, and we denote them with 

I 



Wk, CT (±) = ^ fck )ff +£/) 



±^(e Ka ~e f ) 2 +4\V rT (k)\ 2 D a . 



(50) 
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Following the same arguments employed for the impurity 
we obtain 



n 



^ k,a,t=± 

+N S A(R-l), 



■exp(-/3 w k ,aW] 



where 



o 

fio = -yln 



l + 2exp(-/3e» 



(l + exp(-/3e»)2 



(51) 



(52) 



The same result was obtained by direct analytical inte- 
gration, thus confirming the arguments employed in the 
derivation of Eqs. ( ^l| , pl|) . As in the impurity case, all 
the correlation effects on the thermodynamic potential 
appear in the f2 D , and one expects again a Fermi liquid 
behavior of the quasiparticles in the X-boson approxima- 
tion (cf. the Appendix). 

As in the impurity case we obtain an equation for A 
by minimizing £1 with respect of R: 



A 



S k,cr 

nF (wk, ff (+)) - n F 



■A\V a (k)\ 2 D a 



(53) 



where n F (x) is the Fermi-Dirac distribution. Employing 
|V^(k)| = V as well as a conduction band of constant 
density of states and width W = 2D, we find 



A = — 
D 



D 



n F (wk(+)) - n F (cJk(-)) 



(£k-£/) 2 +4-l/ 2 i? 



(54) 



The quasiparticle properties of a heavy fermion system 
can be described by an effective Hamiltonian, charac- 
terized by two hybridized bands coupled by an effective 
hybridization matrix element Vx3 In the slave boson 
method 



present the same spurious phase transitions shown by the 
slave-boson method, because in all these treatments the 
reduction factor goes to zero at a critical temperature. 
On the other hand, the reduction factor does not vanish 
in the whole range of temperatures for the X-boson ap- 
proach, and the / and c electrons do not decouple. All 
these approaches produce similar results at low temper- 
atures. 



VII. LATTICE RESULTS 



We have employed for the lattice the same system pa- 
rameters used to discuss the impurity, and a similar type 
of behavior is observed. The evolution of the X-boson 
parameters R, and nj are shown in Figure | as a func- 
tion of T, and the / and c electrons do not decouple in 
this treatment because the hole occupation R > in the 
whole range of temperatures, while the corresponding pa- 
rameter r 2 of the usual slave boson treatment vanishes 
at a critical temperature (Kondo temperature Tk)- 
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FIG. 9. Kondo case: Occupation numbers Hf, and parame- 
ters r 2 and R vs. T for both the slave boson and the X-boson 
methods, for the same model parameters used in Figure ^. 



V s _ b = (l-n f )V 2 , 
while in the present approach we have 



V 



X-b 



D a V 2 = (l-n f /2)V 2 



(55) 



(56) 



The hybridization reduction factor is related to the 
effective probability that a c electron jumps into an / 
state, and when U — > oo this transition can only take 
placej4f the / level of the final site is empty. Rice and 
UedaEa argued that the effective hybridization should be 



V 



1 



RU 



1 - n f /2 



V 



(57) 



and variational caleplations were performed by thcmc£l 
and by Fazekas.EJc3 These variational calculations 




FIG. 10. Kondo case: A, and A vs. T for the same model 
parameters used in Figure H. 
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The two approaches yield similar results at low tem- 
peratures. The dependence with temperature of A (slave 
boson) and A (X-boson) is presented in Figure [II], and 
the slave boson A breaks down at the Kondo temperature 
Tk whereas the X-boson A goes continuously to zero, re- 
covering the CHA behavior at high temperatures. 

The evolution of n/, A and r 2 , with the chemical po- 
tential p, is presented in Figure [n] for the slave boson 
treatment. In this case we do not recover the three char- 
acteristics regimes of the PAM: Kondo, intermediate va- 
lence and magnetic, because the formalism breaks down 
when rif — ► 1. The corresponding quantities in the X- 
boson approach: rif, A and R, are also given as a function 
of p in Figure |l2|, and the three characteristic regimes of 
the PAM are clearly shown. There is a plateau in this 
figure when the chemical potential crosses the hybridiza- 
tion gap, and the parameter A is maximum in the Kondo 
region, and goes smoothly to zero together with R when 
\x increases, while fit — > 1. The CHA is recovered in this 
region by the X-boson treatment. 
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FIG. 11. Slave boson: Occupation numbers, A, r 2 , and n/ 
vs. fi for the same model parameters used in Figure H. 
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FIG. 12. X-boson: Occupation numbers, A, R, and nj vs. 
fj, for the same model parameters used in Figure 0. 



The value of the density of states Qf(p), p c {fJ>) as func- 
tion of fi is presented in Figure [L3| for the slave boson 
treatment: the approach breaks down in the Kondo re- 
gion when rif — > 1 and the same quantities are plotted 
in Figure for the X-boson approach. 
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FIG. 13. Slave boson: Pc(p) vs. fj, for the same 

model parameters used in Figure m 
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FIG. 14. X-boson: Pf(fj), p c {p) vs. /i for the same model 
parameters used in Figure H. 



The density of states pf(u>), p c {oj) vs uj are shown 
in Figure [l5] for both the slave boson and X-boson ap- 
proaches in a typical Kondo situation. The density of / 
states at /i is practically the same in the two cases. 
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FIG. 15. Density of states in a typical Kondo situation in 
the two approaches: p{ui) vs u) for the same model parameters 
used in Figure § and T = 0.001. 



VIII. SPECIFIC HEAT 



In this section we apply the theory developed in this 
work to calculate the specific heat of the Kondo insula- 
tors, but for simplicity, we shall not try to fit the ex- 
perimental results of some particular compound. This 
class of materials has been extensively studied in the 
last .decade since its characterization by G. Aeppli and Z. 
FiskEj as Kondo insulators. Some of the compounds that 
we can include in this family are: FeSi; Ce3f3i4Pt3; SmBg, 
and YbBi2. The slave-boson was applied to Kondo 
insulators initially bv-P. RiseboroughcJ and by Carlos 
Sanchez-Castro et al,E3 and we have recently studied the 
Kondo insulator FeSi using the atomic modcl.c^O We 
have been able to adjust simultaneously the static con- 
ductivity, the resistivity and the dynamical conductivity 
to the experimental results, and we obtained a fair agree- 
ment. In this section, we present the X-boson formalism 
as an alternative that gives results close to those obtained 
by the slave boson method in its region of validity, while 
it makes possible to extend the results to the whole range 
of temperatures, without presenting the spurious phase 
transitions appearing in the usual slave-boson method. 

To calculate the specific heat employing ft we first 
show by standard thermodynamic techniques that 



S 



and then that 



OF 
~df 



N.V 



an 
of 



(58) 



C V = T 



dS_ 
~df 



N,V 



-T 



9/i 
df 



' dry 



( ON 



N,V 



(59) 



Assuming a conduction band with constant density of 
states and width W = 2D and the absence of magnetic 
field, we find from Eqs. (51 5^) that 



-T 



d 2 n \ 

dT 2 J 



I3 2 \ - 



= -T 



D 



-T N, 



D 
-D 



d 2 n \ 

dT 2 J 



n,v 



i=± 
d 2 {A(R 



dx u>1(x) rip (tut(x)) [1 — np (ivg(x))] 



1)) 



dT 2 



(60) 



where 



cu ± (x) =\{x + e f )± y( x -s f ) 2 +4\V\ 2 D a 



(61) 



and 



-T 



o 2 n \ 

dT 2 J 



= -2N S k B p 2 ^ exp(/?e» 
[3 + 2exp( / 9e / )] 



[exp((3e f ) + 2} 2 [exp(/3e» + 1]' 



(62) 



We compare the specific heat obtained by the slave 
boson vs. X-boson methods employing in the two cases 
the same parameters, corresponding to a typical Kondo 
insulator situation with the chemical potential inside the 
gap. In Figure 16 we present the corresponding density 
of states Pf{oj) vs. uj, for the following parameters: Et = 
0.3; V = 0.35; W = 2.0; T = 0.001 and /i = 0.5. 
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FIG. 16. Density of states in a typical Kondo insulator 
situation in the two approaches: Pf(u>) vs. ui with the fol- 
lowing parameters: E t = 0.3; W = 2.0; V = 0.35; fi = 0.5; 
T = 0.001. 

In Figure [Tt] we present C v vs. T, employing the 
same parameters of Figure [T|. We have calculated 
T (dfi/dT) N v (dN/dT)^ v numerically, and its contri- 
bution to Cy is negligible for these parameters. The cal- 
culation above the Kondo temperature Tk in the slave 
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boson case was performed for the phase of uncoupled / 
and c electrons, and the specific heat presents a discon- 
tinuity at the transition Tk, but the curve obtained by 
the X-boson method is continuous. 
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FIG. 17. Kondo insulator: Specific heat of a typical Kondo 
insulator for the same model parameters used in Figure |l6| 
but as function of temperature. 



IX. CONCLUSIONS 

The slave-boson formalism! has .he.en extensively used 
in the mean-field approximatior£3o to study strongly 
correlated electron systems, because it captures the es- 
sential physics of the Kondo effect at low temperatures 
and its implementation in the Kondo limit is very simple. 
One drawback of this method is that above a tempera- 



ture Tk (cf. Section IV) or when fi » Ef, it develops 
a spurious second order phase transition into a phase 
with decoupled conduction and localized electrons. In 
the present work we present an approach that was in- 
spired in some of the features of the slave-boson method 
in the mean-field approximation and gives similar result 
at low temperatures, but it does not have the spurious 
phase transitions of that method at higher temperatures 
or chemical potentials, where it behaves like the CHA 
solution. We have named "X-boson" this new technique 
because it keeps the X operators!] even at the final stages 
of the calculation, while two hybridized but uncorrelated 
bands appear in the slave-boson method at that stage. 
What we employ instead is an approximation to the exact 
cumulant expansion of the PAM in the limit x)f U — > oo, 
which is called CHA (chain approximation) .Era The CHA 
corresponds to the most general set of diagrams with only 
second order cumulantaj, and has interesting properties 
like being <!>— derivable. B We have developed the X-boson 
method both for a single local impurity and for a lat- 
tice. The correlation appears in this approximation only 
through the presence of a constant Do a at critical places 
of the corresponding GF (cf . Eqs . (|l| - ||, || - [50|) ) , but 
it has profound effect on the system description, forcing 



le, the C p oiFeSi measured 



the local electron occupation to be nf < 1, and elimi- 
nating the spurious phase transitions of the slave-boson 
method. 

From the computational point of view the X-boson 
method is equivalent to the slave-boson technique, and 
the possibility of using it for all temperatures or chemi- 
cal potentials fi makes it useful as the starting point to 
study several interesting problems: competition between 
Kondo eifjepi. and RKKY interaction in heavy fermion 
sy stems, Ejlj impurity bands in Kondo insulators, til 
non Fermi- liquid behavior in disordered systems ,E3 and 
Kondo effect in quantum dots.BE-a To give an example 
we have considered a Kondo insulator and calculated its 
specific heat: the slave boson method shows a discontinu- 
ity at Tk, while it is continuous in the X-boson method. 
Although it is C p that is usually measured, the difference 
C p — C v is usually small in liquids and solids, and shows a 
dependence with T similar to that obtained with the X- 
boson method (see for exa: 
by Jacarino and coworker! 

In conclusion, we have presented an approach that 
gives essentially the same results as Coleman's slave- 
boson treatment at low temperatures, but without the 
spurious phases transitions at intermediate temperatures 
T or large values of /i. The X-boson results approach 
those of the cumulant expansion in the chain approxima- 
tion for large values of T or fi. 
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APPENDIX I: FERMI LIQUID BEHAVIOR OF 
THE X-BOSON METHOD 



We show here that the PAM has Fermi liquid proper- 
ties in the X boson description, and to discuss the prob- 
lem we shall first consider the / and c electrons self- 
energies S/(k, z) and S c (k, z), that are related to the 
exact GF by the equations 



(z -e), ff -£/ (k, z)) 



(z - e^.a - S c (k, z)) ' 



(1) 



(2) 



where we employ the z with the y, already subtracted, 
(cf. Eqs.( ^0| , pl] )). Employing an effective cumulantc3'EZl 
M^J{z) it is possible to write the exact GF 



G f d(z) 



Mf ff 



1 2M 



l-\V(k)\*G°,(z)Me"(z) 



(3) 
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and 



It then follows that 



l-|nk)| 2 G° kCT (z)M 2 e ^(z) 



C.ktT v 

-1 



z- £kff + |y(k)pM^(z) : 



where 



(4) 



(5) 



is the unperturbed GF of the c electrons. The exact self 
energies are then 



£,(k,z) 



and 



|y(k)| 2 G° kCT (z) 



S c (k,z) = -|y(k)| 2 M 2 e ^(z), 



(6) 



(7) 



and to recover the CHA self energies it is sufficient to 
replaceB the effective cumulant Ml -ISz) by the contri- 
bution of its simplest diagram, namely the unperturbed 
GF of the / electrons: 



AL 



eff 
2,ker 



{z)^G° f)0<T {z) = -D 0a /(z-e tF ). 



(8) 



The essential property of a Fermi liquid is that at low 
temperatures it has a behavior close to that of a system 
of free Fermions. It is described by a system of quasi- 
particles that replace the elementary excitations of the 
free system, but they have a finite life-time caused by the 
interactions that are absent in the free system. Instead 
of the elementary particle energies we have to analyze 
the poles of the relevant GF, given in our problem by the 
solutions Zko- of the equation 



e^ + \V(k)f M% a (z) = 



(9) 



The Re [z-^ a ] corresponds to the energies of the elemen- 
tary excitations, while their imaginary parts give their 
decay properties, and the Fermi surface is given by the 
set of kp such that Re [zka] — 0, because we already have 
subtracted the /i from the z. The M%-J^{z) in the nu- 
merator of Eq. does not introduce new poles, because 
G(/ a (z) is finite at those poles. 
/From Eq. (Bj) we obtain 



|V(k)|* ' 
and replacing in Eq. (^|) we find 

S/(k, Z k a) = Z kcr - £/,<x, 

while from Eqs. (|]||) we obtain 

E c (k, Zkcr) = ^kCT — £k(7- 



(10) 

(11) 



Im [ E c (k, z kCT )] = Im [ £/(k, z kcr )] = Im [z krT ] . 



(12) 



An essential property of the Fermi liquid is that the 
quasi-particles on the Fermi surface have an infinite life- 
time when T — > 0, and this property follows from Eq. 
( |l2| ) for both the / and c electrons when the poles on the 
Fermi surface are real, because the corresponding self- 
energies are then also real. 

/.From Eq. J5C| ) we conclude that the GF poles of 
the PAM in the CHA are all real, and the same prop- 
erty holds for the X-boson solution. These two approx- 
imations therefore describe the PAM as a Fermi liquid, 
because they satisfy the condition above. Because the 
infinite lifetime is valid for all k and for all T, the quasi- 
particles in these two approximations behave more like 
elementary excitations than like the usual quasi-particles, 
that have a finite lifetime outside the Fermi surface, even 
at T -> 0. 
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